Theory of AC Anomalous Hall Conductivity in d-electron systems 
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To elucidate the intrinsic nature of anomalous Hall effect (AffE) in d-electron systems, we study 
the AC anomalous Hall conductivity (AHC) in a tight-binding model with (d^z, dj;2)-orbitals. We 
drive an analytical expression for the AC AHC axy{uj), which is valid for finite quasiparticle damping 
rate 7=ft/2r, and find that the AC AHC is strongly dependent on 7. When 7 — +0, the AC AHC 
shows a spiky peak at finite energy A that originates from the interband particle-hole excitation, 
where A represents the minimum band-splitting measured from the Fermi level. In contrast, we 
find that this spiky peak is quickly suppressed when 7 is finite. By using a realistic value of 7(cij) 
at oj = A/2 in d-electron systems, the spiky peak is considerably suppressed. In the present model, 
the obtained results also represent the AC spin Hall conductivity in a paramagnetic state. 

PACS numbers: 72.10.-d, 72.80.Ga, 72.25.Ba, 72.25.-b 



I. INTRODUCTION 

In usual metals, the ordinary Hall effect due to Lorentz 
force is widely observed. In this case, the Hall resistivity 
is proportional to the applied magnetic field. In ferro- 
magnets, in addition, the Hall resistivity that is propor- 
tional to the magnetization is observed, which is called 
the anomalous Hall effect (AHE). Therefore, the conven- 
tional expression for the Hall resistivity, p^, is given by 
PH = R^^B + iTTRjfM, where R% and R% are the ordi- 
nary and anomalous Hall coefficients, B is the magnetic 
field, and M is the magnetization. The mechanism of the 
AHE has been intensively studied for a long time. 

The study by Karplus and Luttinger (KL) 1] in 1954 
was the first theoretical approach to the AHE, and it 
was refined by Luttinger [3] in 1958. They pointed 
out that the anomalous Hall conductivity (AHC) <J^y{= 
R'^M / p^) is finite and dissipation-less {cr%y oc p^) when 
Af ^ in multiband systems with the spin-orbit interac- 
tion (SOI). This KL-term is called the "intrinsic AHE" 
because it is due to the interband particle-hole excita- 
tion, and it exists even in systems without impurities. 
On the other hand, alternative mechanism due to impu- 
rities was proposed as the "extrinsic AHE". Smit have 
shown that extrinsic AHC due to the impurity skew scat- 
tering follows crjj; cx p^^ Later, Berger found that 
another extrinsic mechanism, the side jump mechanism, 
gives a'^Jy oc Recently, AHE due to skew scat- 

tering and side jump mechanism has been studied based 
on linear response theory and semiclassical Boltzmann 
equation approach d, H, 0] ■ 

After KL, intrinsic AHE has been studied by several 
specific theoretical models d, 0, i, i, [13, El, [ll [H, Hi, 
m, [3. Recently, AHE in the Rashba 2D electron sys- 
tems has been intensively studied [13, [H, [H, H [SI, [Ijl . 
In refs. [20, iH, [H, they have reported that the AHC 
vanishes in the Rashba model due to the cancellation by 
the current vertex correction (CVC) due to impurities 
unless the lifetime is spin-dependent. In graphene sys- 
tem, the appearance of large quantum spin Hall effect 
has been predicted when the Fermi level lies inside the 



gap (23|,|24|,[25|. 

In general, intrinsic AHC is composed of the "Fermi 
surface term" and the "Fermi sea term" [l^]. Recently, 
refs. [13 , list have reported that KL's AHC, which is a 
part of the "Fermi sea term" , is expressed in terms of 
the "Berry curvature" in case of 7 = 4-0, where 7 is 
the quasiparticle damping rate. In several simple models 
!l2„ JJi, J^,i23j], the Berry curvature term gives the correct 
AHC since another Fermi sea term exactly cancels with 
the Fermi surface term. However, it presents erroneous 
result in the high resistive regime since the cancellation 
of other terms becomes imperfect, [l^, [27l [28j. 

AC transport phenomena have been attracting much 
interest for a long time since they can provide detailed 
and decisive information on the magnetic properties and 
the electronic structure of magnetic materials. For exam- 
ple, the AC AHE have been studied as the well-known 
magneto-optical effect (MOE) . In transition metal ferro- 
magnets, theoretical studies of MOE had been reported 
in refs. [H, [lO, |4lj. Therein, the overall behavior of 
the experimental AC AHCs at high freqencies are repro- 
duced in many transition metals. Moreover, in high-Tf 
superconductors (HTSCs), AC Hall effects have been in- 
tensively studied by Drew's group. They have found that 
the Hall cocfBcicnt Rh{^) in HTSCs shows a prominent 
w-dependence |2i, |30|, |3l|, [3^, Q IM |35| , which can not 
be reproduced by the extended Drude form. This anoma- 
lous AC transport phenomena can be explained by the 
fluctuation-exchange (FLEX) approximation by consid- 
ering the CVC [Sg [stI. . These studies on AC trans- 
port phenomena showed the significance of the strong 
antiferromagnetic fluctuation in HTSCs. In similar, AC 
AHE at low frequencies will be a fruitful study to under- 
stand the intrinsic nature of AHE. 

Recently, AC AHEs in SrRuOs and Fe were calculated 
based on the LDA band calculations [ij, [l^ . According 
to the ab initio calculation [i^], the AC spin Hall effect 
(SHE) for p-type semicondutors was also studied. In ref. 
|14S|, they have predicted that AC AHC has a sharp and 
spiky peak at low frequencies that originates from the 
interband transition under the assumption that 7 — +0. 
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However, in metallic systems, the damping rate 7(1^) 
should be finite for w ^ even at zero temperature be- 
cause of inelastic scattering: In a Fermi Liquid, the quasi- 
particle damping rate is given as 7(w) oc {(ttT)^ 
where T represents temperature. Therefore, a reliable 
calculation on the AC AHC for finite 7 is highly required. 

The aim of this paper is to obtain the reliable expres- 
sion for the AC AHC for finite 7 based on the {dxz, dyz)- 
orbital tight-binding model. For this purpose, we derive 
an analytical expression for the AC AHC including both 
the Fermi surface and Fermi sea terms based on the lin- 
ear response theory. This expression is valid for finite 7 
if the CVC is negligible. Using this expression, the AC 
AHC at low frequencies, where uj < 4000K, is studied in 
detail. We find that the intrinsic AC AHC shows a non 
Drude-like behavior. In the case of 7 = -1-0, it has a sharp 
and spiky peak a,t u; — A. Here, A is the minimum band- 
splitting measured from the Fermi level, and A > lOOOK 
in usual transition metals. However, this spiky peak is 
easily suppressed by finite 7. By using a realistic value 
of 7(w) at w = A/2 in usual metals, the spiky peak is al- 
most smeared out. We also find that the overall behavior 
of the AC AHC is reproduced well by the Fermi surface 
term, whereas the Berry curvature term reproduces the 
correct AC AHC only when 7 ^ A. Since this condition 
is not satisfied in usual metals for uj ^ A, the Berry cur- 
vature term gives erroneous AC AHC in the real metallic 
systems. 

Now, we explain that AC AHE at low frequencies can 
provide important information on the mechanism of the 
AHE. For a long time, the origin of AHE (intrinsic or 
extrinsic) in real systems has been unsettled problem. 
In clean heavy fermion systems, a^y is independent of 
p (i.e. Rh oc p^) sufficiently below the coherent tem- 
perature Tn, whereas a^y oc {Rh oc p°) above Tq 
[il I45I lia |47| . This fact indicates that the intrinsic 
AHE is dominant in clean samples [1]. Also, recent ex- 
periments for several transition metal complexes have 
reported that tr" is independent of p in the low resis- 
tive regime psl |49| , whereas a^y decreases in proportion 
to p~^ in the high resistive regime [49]. Based on the 
intirinsic AHE in (d^z, dyz)-OThital tight-binding model, 
Kontani et al. have explained this experimental re- 
sult in transition metal ferromagnets. However, in DC 
AHE, it is not easy to distinguish experimentally which 
mechanism (intrinsic or extrinsic) is dominant since the 
introduction of randomness, by which both the resistivity 
p and the skew scattering increase, makes the analysis of 
experimental results difficult. In contrast, AC AHE may 
be useful to solve this problem by measuring disorder 
free samples: Since the AHC due to the skew-scattering 
mechanism is proportional to 7"^ if elastic scattering is 
dominant, the AC AHC due to the this mechanism shows 
a sensitive w-dependence like the Drude-type behavior: 
cr^y{uj) oc (7 — jcj)"^ for small cu. On the other hand, 
AC AHC due to the intrinsic mechanism is independent 
of lli for (jj ^ A, as shown in this paper. Therefore, AC 
AHE will be useful to resolve the controversy over the 



origin of the AHE without necessity to introduce disor- 
ders. 

Finally, we comment on the AC SHE: SHE is the phe- 
nomenon that an applied electric field induces a spin cur- 
rent in a transverse direction. Recently, refs. [13, [H, l43j 
have found that the huge SHE is ubiquitous in multi- 
orbital d-electron systems. The origin of the huge SHE 
is the "effective Aharonov-Bohm (AB) phase" induced 
by the atomic SOI with the aid of inter-orbital hopping 
inegrals. Since s^-spin current operator j| is given by 
— — I (ja;| — jxi) in the present model, the relation 
a^y = —^(Txy holds. Here, jx^Uxi) is the charge current 
operator for f (J,)-spin, and a^y = j^/Ey. Therefore, in- 
teresting w-dependence of AC AHC derived in the present 
study is also expected to be realized in AC spin Hall con- 
ductivity (SHC). 



II. MODEL AND HAMILTONIAN 

In this paper, we study a square lattice tight-binding 
model with dxz- and dyz- orbitals, which is a simplified 
model of a famous triplet superconductor in Sr2Ru04 
[sot . The {dxz,dyz)-OThital tight-binding model is one of 
the simplest models for studying the AHE in transition 
metal ferromagnets [T6|. Using this model, we derive 
explicit expressions for AC AHC that is valid for finite 
quasiparticle damping rate. 

To realize the AHE in ferromagnets, the SOI is indis- 
pensible. In a perfect ferromagnetic metal with A^Hz, 
the atomic SOI AZ • J is replaced with —X{A'I/p,b)Iz, 
where A is the coupling constant. The intrinsic AHE 
is caused by the interband transition of quasiparticles 
due to the off-diagonal elements of Iz [H. In the d- 
electron systems, the matrix element of Iz is finite only 
for {yz\lz\zx) = -{zx\lz\yz) = i and {xyllzlx"^ - y^) = 
— {x^ — y'^\lz\xy) = 2i. Note that \zx) and \yz) are 
given by the linear combination of = ±1, and \xy) and 
\x'^ — y"^) are given by the linear combination of Iz = ±2. 
In pure transition metals, since the energy splitting be- 
tween t2g and Cg orbitals is smaller than the bandwidth, 
the interband transition of the quasiparticle between dxy- 
orbital (in t2g) and dx^-y^-orhital (in e^) should be taken 
into account correctly. In fact, we have found that the 
SHE is mainly caused by {dxy , dx^ _y2)-OThital in pure 
transition metals ^28j . On the other hand, in ruthenates, 
since the energy splitting between t2g and Cg orbitals are 
large, AHE is maily caused by {dxz, (ij,z)-orbital. In fact, 
{dxz, dyz)-OThital tight-binding model can explain impor- 
tant experimental fact as reported in ref. [4^. We note 
that the mechanisms of large AHEs arises from {dxz, dyz)- 
orbitals and {dxy, (ia;2_j,2)-orbitals are the same: The "ef- 
fective AB phase" factor of conduction electrons due to 
d-atomic angular momentum with the aid of the atomic 
SOI and the inter-orbital hoppings [2l[2i,[4l|. 

Here, we represent the creation operator of an elec- 
tron on dxz-{dyz-) orbital as c^^(c^^). The Hamiltonian 
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without SOI is given by i?" = c^^/i^Cfc, where [Tf 



(1) 



and cl 



-2i cos fcj^ and 



^ = 4i' sinfca; sinfcy, where —t and ±t' are the hopping 
integrals between nearest-neighbors and nextnearest- 
neighbors, respectively. They are shown in Fig. [T] 

In the present model, the velocity matrix — 
dh^/dk^ {fj, = x,y) is given by 




FIG. 1; (a) Hopping integrals between the same orbitals. 
(b) Hopping integrals between the different orbitals, which 
is given by the nearest neighbor hopping. Note that the sign 
of the interorbital hopping integrals change by 7r/2 rotaion. 
This fact gives rise to the anomalous velocity. 



2t sin kx 
At' cos kx sin ky 



At' sin kx cos ky 



At' cos kx sin ky 


At' sin kx cos ky 
2t sin fc„ 



(2) 
(3) 



We should stress that the off-diagonal elements of Vx , 
are odd- functions of ky. In the same way, 



Vy^ = v"^^ are odd-functions of kx- They are called the 
"anomalos velocity". In later sections, we will see that 
AHC axy is proportional to \{v%^Vyy). This implies that 
sizable AHC is caused by anomalous velocity with the aid 
of atomic SOI. Consequently, atomic d-orbitals degree of 
freedom gives rise to the huge AHC in transition metal 
ferromagnets. 

Next, the atomic SOI, which is indispensible for AHE, 
is given by Ti/'^ = Y^ t. c]M^Ck- Here, in the present 
bases is given by [If 



-sgn(s^)Afj^, 



(4) 



where Ty represents the Pauli matrix for the orbital space. 
Hereafter, we put = 1 for the simplicity of calculation. 
The Green function in the presence of atomic SOI is 

given by Gk{'^) — + ^ h'^^ , which is ex- 

pressed as follows in the present model [l^ : 



r>i, A — f ^x^i^) Gxy{^) 

^^'^^ - 1, Gyx{^) Gyyiu;) 



d{uj) 



Oik 



(5) 
(6) 



where ak = Cfc^ ~ i^sgn{sz) and d^iv) = 
/I — — |Q!fep, which is expressed as 

d{uj) ={lu + /i - E+){uj + E^), 



(7) 
(8) 



Here, represents the quasiparticle dispersion. Figures 
[2] (a) and (b) show the Fermi surface and the band struc- 
ture obtained in the present model for {t,t') = (1,0.1), 
which is realized in Sr2Ru04 [lH. The electron density 
per spin, n, is set as 0.4. The energy splitting rep- 
resents the minimum band-splitting {\E+ - E^\) mea- 
sured from the Fermi surface of E'^-band. In Fig. [21 
= \E'^ — Ej^ \ at k'^, and k* represents the position 
of the minimum band-splitting A = minjA"*", A~}. 

Here, we consider the quasiparticle damping rate r(ci;). 
Microscopically, it is given by the imaginary part of the 
self-energy, For simplicity, we assume that 

is diagonal with respect to orbitals, and is independent of 
the momentum: Tap{uj) — ^{uj)5ai3, where a and (3 are 
orbital indices. This assumption is justified in the case 
of the elastic scattering due to local impurities, since the 
local Green function gap = Gap is small for a ^ (3 
in the case of A <C Wband, where Wband is the bandwidth 
[l^. This fact is also justified in the case of inelastic 
scattering due to on-site Coulomb interaction in the dy- 
namical mean field approximation (DMFA), where the 
self-energy is composed of local Green functions. Using 
eq. ([6]), the retarded and advanced Green functions are 
given by 



G^p,{uj) = Gap'{uJ + ijiu;)), 

GaP'i^) = GaP'{uJ - il{uj)), 



(9) 
(10) 



respectively. Hereafter, we assume that the renormaliza- 
tion factor Zk={l~ dt^/dcuj = 1 since it cancels out 
in the final expression of the AHE [l^ . 



III. 



AC HALL CONDUCTIVITY 



In this section, we derive an analytical expressions for 
the AC Hall conductivity axy{uj) and the AC longitu- 
dinal conductivity (Txxi^) based on the linear response 
theory by dropping all the CVC. The other approach to 
the AHE, which is the semiclassical Boltzmann equation 
approach, is also useful d, 0, [13: In ref. [1], the equiv- 
alence of these two methods has been shown in the 2D 
Dirac-band graphene system. 

Until section |Vl we assume that the elastic scatter- 
ing due to the impurity potential is dominant over the 
inelastic scattering due to electron-electron interaction. 
In this case, the quasiparticle damping rate is given by 



-, where N{lo) is the density of 



1 - t:IN{uj) 

states per orbital, riimp is the impurity concentration. 
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FIG. 2: (a) Fermi surface for n = 0.4. Here we put {t, t') — 
(1,0.1), which corresponds to Sr2Ru04. Here, F — (0,0), 
M = (tt, tt), and X = (0, tt). (b) Band structure in the present 
model. Here, we set n=0.4, A=0.2 and t/t' — 0.1. 



and / is the impurity potential. Since the oj-dependence 
of N{uj) is small in the present model, we omit the uj- 
dependence of jimp- In section |Vl we study the AC 
Hall conductivity in systems where inelastic scattering 
is dominant. In this case, we have to take account of uj 
dependence of the damping rate. 

Here, we comment on the CVC due to the local im- 
purity potential. In the Born approximation, the lowest 

order CVC is given by AJ^ = rumpl'^ ^^G'^Jf.G^ ■ 



In the d-orbital tight-binding models, G has an even par- 
ity with respect to fc ^ — fe. Therefore, AJ^ = since 

In 



{d/dkf,)G = GJ^G is an odd function [16|, |28|, [42 
contrast, the CVC plays an essential role in the Rashba 
models [S, HH, [H • 



According to linear response theory, the AC Hall con- 




FIG. 3: Diagrammatic expression for axy{to„) 

ductivity is given by 

a,A^) = ^[K^U^)-K^'M], (11) 

where K'I}^{lo) is the retarded correlation function, which 
is given by the analytic continuations of the following 
thermal Green function: 



k,k' ,a^a' ,13,13' 



vZ^^vt^Qkk'ioj^), (12) 



where 



Vk.t 



dkfj.' 



(13) 



and 



Qkk'M = J dTe'"""(r, (aL'(^)afca(T)4';3'(0)afc'^(0)))- 

(14) 

If we drop the CVC, Qk,k' is given by 

Qkk'{^n) = —6kk'P ^ ^ G/3q'((jJ„JG'q'/3(w„ + CJ„J, 
ni 

(15) 

The diagrammatic expression of (J^^{LUn) is shown in Fig. 

m 



Performing the analytic continuation carefully, the expression for the AC Hall conductivity axy{uj) is given by 

k.a' ,a.l3' ,p 

X [{fie + u;/2) - fie - c./2)} Gf^, (e + ^/2)G^„,(£ - c./2) 

+/(£ + Lo/2)G^^, ie + Lo/2)Gf^, (e - c./2) - fie - u:l2)G%. (e + c./2)G^„, (e - t^/2)] (16) 
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where f{e) is the Fermi distribution fuction. Here, we divided axy{i-^) into two terms, where the first term in the 
square bracket corresponds to the Fermi surface term (/), and the second and third terms correspond to the Fermi sea 
term (//). Since / term consist of G^G^ , whereas // term consist of G^G^ and G^G"^, the division of axy{uj) into 
these two terms is unique. Hereafter, we drop the factor e'^/h in a^y{uj) to simphfy expressions. Note that K^y{Q) = 
since it is proportional to G^G^ - G^G^ [l6j . 

Now, we first take the summation over a' ,a, (3' , (3 in eq. p6p . In the present model, the terms {a' , a, (3' , (3) = 
{x, X, x , y) , {x, X, y, x), (x, y, x, x) and (y, x, x, x) remain finite after fc-summation, just like the calculation of DC AHC 
in ref. [l6|. Considering the square lattice symmetry of the present model, we obtain the following expression: 



k 

X {f{e + Lu/2)- f{e-uj/2)} 
+fie-oj/2) 



oj + 2ij 



d^{e + uj/2)d^{e-cu/2) 



d^e + uj/2)dR{e-uj/2) 



fie + u/2) 



d^{e + uj/2)d^{e~uj/2) 



(17) 



This integration by e can be calculated analytically, and the final result for the AC Hall conductivity axy{uj) at T—0 
is given by 



1 



El~E-~u;^2i^ 
1 

Ek -Ek +UJ + 2ij 



In 



In 



fi + iv - + Z7 

l^-Et - il 
H + LU - + i-f 



In 



In 



^-uj- E+ -ij 



■ 17 



Et-El 



Et-El-r^ 



{In - E+ + Z7) + In (/i - E^ - i^) } 



Et-E^^-u: 



-iX 



1 



xy 

' Et-El 



Et - E- 



Et-E^-u: 
- i X [cq.®]. 



- {In + cj - E^ + i-i) + In (/i - - - ^7) } 
{In + - £;+ + i-f) +\n{^l-uJ-E- - i-f)] 



(18) 



(19) 



(20) 



(21) 



If we put a; = 0, the obtained expression reproduces the DC AHC given in ref. [lg|. In the above equations, we 
divided the Fermi sea term into two terms, which are aj^{Lo) and cr^^^(w): cr^^'' is called the Berry curvature term 
[1^ IT3I [4^ [5^. Here, we analyze eqs. ([T9|) . ([20]) and 0A\ in the clean limit n,;„ip ^0 (7^0), and show that the 
real part of aly^{uj) in eq. corresponds to the expression which are frequently used to calculate (Txy{i^)- In the 
band-diagonal representation {E^; a = ±), we can show that the off diagonal velocity is given by 



XX xy I „,xy yy 
„,aa aa _ _ ■ ^ ''x ""y ^ "^x "y 

^ ^ ~ EZ-Et 



(terms which vanish with fc-summation.). (22) 
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Since Im [ln(2:) ± ij] = ^tt6{x) for 7 — * +0, eq. ([20]) is rewritten as follows in the band-diagonal representation: 

a,, (a;) = }^ (23) 

where m, n are the band indices, and an infinitesimal imaginary part i5 is included to maintain the analytic property 
of the Hall conductivity. The obtained expression of <Tly''{Lu) in eq. (|23p corresponds to the expression which are used 
to calculate axy (w) in refs. [s^, HO, Sil, [42, [S^ . 

In literatures, eq. ([^5]) has been frequently used since aly{uj) + cr^^" = for 7 = +0. However, this condition for 7 
will not be satisfied in the real metallic systems since 7 increases with uj due to inelastic scattering. For this reason, 
axy{oj) ^ al.^y{uj) in usual metals. Therefore, for reliable calculation, we have to analyze I, I la, and lib terms on the 
same footing. 

Here, we comment on the physical meaning of the three terms: a^y and c^y" are finite only when the Fermi surface 
exists, whereas cr^^'' is finite even if the Fermi surface is absent [l^. a^l'' term is the origin of the quantum Hall effect 
and the quantum spin Hall effect [H, [13, [H, [H, H [SS^ . 

According to eq. (|lip . the expression for the AC longitudinal conductivity axxii-^) is given by 



27r " " 



k.a' ,a,(i' ,(i 

X { /(e - u,mG%, [e + C./2) {G% (e - u;/2) - Gj^, (e ~ cj/2)) 
+ f{e + oj/2) {G^p, (£ + uj/2) - Gtp, {e + c./2)) G^„, (e - a;/2) 

-/(£) {G^p.{e)Gf,Ae) - G%.{e)G^^,{s))] . (24) 

In the same way as axy{^), we can calculate axx{^) analytically. By dropping the terms that vanish after fc-sumniation, 
axxiyj) is given by, 

'^-(^) = [{vT?A{u^) + 4«r<^i?(a;) + 2{v:yfC{u)\ , (25) 



where 



and 



' 27ra; 

k 



A{u:) = /12 - 2(el - ^7)/ll " {(^/2)^ - (Cl " ^if] ho 

[h, 2{a + ^l)hl {{i^/2f ia + n)'} /30] 

- [I22 - 2^121 {(C./2 + (^D'} ^20] 

- 2ilm [K2 ~ 2(e^ - ij)Ki + - i7)'i^o] , (26) 

B{co) = ^r^Ai - a'ia - n)ho - i^i^hi - + *7)/3o] - r^'hi - ^i^ekU 

-2^C^lm[K,-{a-^J)Ko], (27) 
C(c.) = /12 - iCk +a- 2^7)/ii + {(O' - + (e^ - Z7)(^^ - Z7) - (c./2)2} 

- [h2 (S + + 2.7)/3i + - + (e^ + .7)(e^ + ^7) - (c./2)2} /30] 

- [122 iCk + a)i2i + { iCk^r - + m (-'/2 + *7)'} /20] 

- 2zlm [K2 + a- 2.7)^1 + {{ek'f - A^ + (e^ - - *7)} i^o] , (28) 

hn = / cJe-TBT \ ,p/ = / C^E" 



— 00 



' dR{e + uj/2)dR{e-u/2y 7_^/2 rf^(£ + c^/2)d^(£ - cc;/2) ' 



0/2 1-0 



^'-=1 '^'d^{s + ul2)d^{e - U./2) ' ^" - i_ '^^d«(£)d«(e) ("^O'^'^)- 



00 



These integrals can be performed analytically. As a result, we obtain the expression for the AC longitudinal conduc- 
tivity axxii^) at T—0. We can verify analytically that this expression reproduces the longitudinal conductivity axx at 
uj = 0, which is given in ref. (l6j . 



IV. NUMERICAL STUDY feet ferromagnetic state where = n and 

In this section, we perform the numerical study for 
both axy{ijj) and axx{'^) at T = 0, assuming a per- 
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In this case, niz = I^bti. Hereafter, we put = 1. 
Also, we put the couphng constant of SOI as A = 0.2. 
It corresponds to 800K if we assume t=4000K, which 
is a realistic value in ruthenates. The main purpose 
of this section is to elucidate the frequency (lu) depen- 
dence and the damping rate (7) dependence of the AC 
Hall conductivity. We perform the fc-summations in eq. 
psp for (Jxy and in eq. (|25p for axx numerically, di- 
viding the Brillouin zone into 5000 x 5000 meshes. By 
using the obtained results for a^yiuj) and a^xi^^), we 
also present the uj- and 7-dependences of the Hall co- 
efficient Rh{^^) — <^xy{^jj)/<y'ix{'^) the Hall angle 
Onii^) = (Jxy[uj)/axx{i^)- 

The unit of conductivity in this section is /ha, where 
h is the plank constant and a is the unit cell length. If 
we assume the length of unit cell a is 4A, then /ha « 
lO^fi-icm-i. 

Figure [3] shows the cj-dependence of Re axy{oj) for 7 = 
0.001,0.1, and 1. When 7 = 0.001, Re (Jxy{oj) has a 
sharp peak at w A as shown in Fig. [4] (a). In this case. 
Re ax^y{Lo) (Berry curvature term) reproduces the total 
AHC Re axy{uj) well. On the other hand, we see from 
Fig. |4] (b) and (c) that the spiky peak of Re axy{(jj) is 
significantly suppressed when 7 is large. This fact is well 
reproduced by the Fermi surface term CT^y(a;), whereas a 
sharp peak remains in the Berry curvature term (t^^^(w) 
against large 7. Therefore, al.^y{uj) gives incorrect reslut 
of AHC when 7 is finite. We note that another Fermi sea 
term aly°'{uj) is not shown in these figures since two Fermi 
sea terms satisfy the relationship aly'^{uj) ~ — cr^^^(w) for 
any value of 7. As a result, the Fermi sea term is small 
in magnitude, and axy{uj) gives a main contribution to 
axy{(^)- Thus, (Txy{i^) repoduces the total AC AHC well 
for wide range of to: 



axy{(^) ~ (T^j,(ijj) (Fermi surface term). 



(29) 



This is one of the most important results in the present 
paper. Recently, the LDA calculations of AHC and SHC 
in real systems were performed by considering finite 7 
[5^. [57I. \5§\ . However, the effect of 7 is underestimated 
in their calculations since only cr^^'' is calculated. 

We have checked the reliability of the numerical study 
in two ways: First, we verified that the AC conductivities 
(aU axx{uj),a-xyiuj),aiy{uj), ail^iuj) and cr^^''(w)) repro- 
duces the DC conductivities (w = 0) derived in ref. \W( . 
Second, we calculated the sum rule for Oxyito) numeri- 
cally. The sum rule for axy (w) is given by 



I do; Re axy{'jj) — 0. 
Jo 



(30) 



Equation ([50)) is easily recognized from facts that 
Gxyiyi) ~ as — > 00, and it is ana- 

lytic in the upper-half plane of the complex w-space. 
Here, we performed the numerical w-integration of 
r = divRe Oxyioj)/ /J^dw|Re cra;j,(w)|, where we put 
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FIG. 4: cj-depenence of Re axy{(jj), Re aiy{Lo) and Re aiy''{Lo) 
for 7 = 0.001 in (a), for 7 — 0.1 in (b), and for 7 = 1 in (c). In 
(b) and (c), a sharp peak at ~ A in Re ali/'{uj) is erroneous. 



X=1QQ. It should vanish identically when X = 00 ac- 
cording to the sum rule. We have verified that r ^ 10~^, 
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which suggests the high rehabiHty of the present numer- 
ical study. 



A=0.417 



between DC AHC and Re (Jxy 



UJ 



It is a well known 




A=0.417 



'pO.OOl 
'pO.Ol 



-pO.05 

y=0,\ 

•l y=0.5 




FIG. 5; cj-dependence of axy{uj) for various 7. The real part 
of axy{i-o) is shown in (a), and the imaginary part in (b). The 
spiky peak in Re axy{i^) at ~ A is quickly suppressed by 7. 

Hereafter, we show the numerical results of AC AHC 
in more detail. We first discuss the 7-dependence of the 
AC Hall conductivity a^yioj)- Re a^yiuj) for various val- 
ues of 7 is shown in Fig. [5] (a). When 7 is very small, 
Re iJxy{uj) has a sharp peak at finite energy A. After 
reaching the peak at w ~ A, Re axy{uj) decreases drasti- 
cally and changes its sign. According to eqs. (fT8|) . ((20|) 
and (|2ip . the main contribution for Re a^yi^) comes from 
area near k* in Fig. O When 7 becomes large, however. 
Re axy{id) becomes almost constant for < < A, and 
the peak at w ~ A vanishes. In Fig. [5] (b), we show the 
w-dependence of Im axy{u>) for various 7. For 7 = 0.001 
and 0.01, Im axy{io) starts to increase drastically around 
a; = A, and it takes a maximum value at w ~ 0.5 which 
is slightly larger than A. We see that this peak is sup- 
pressed as 7 increases. 

Now, we discuss the difference of the 7-dependences 



property that DC AHC cTxy is independent of 7 in the low 
resistive regime where 7 ^ A [HH, This property 
can be recognized in Fig. O (a) at w = 0. In contrast, 
at finite frequencies. Re axy{u!) for 7 = 0.001,0.01, and 
0.05 in Fig. \^{a) behaves quite different from each other, 
especially around w ~ A. Thus, the 7-dependence of 
the AC AHC is much more sensitive to the value of 7 
compared with the DC AHC. 



A=0.417 




FIG. 6: (^-dependence of Re <Jxx{io) for 7 

0.001,0.01,0.05,0.1 and 0.5. 



We also discuss the w-dependence of Re axx (t^) . Figure 
[S]shows that Re axx has the Drude peak at = 0. It also 
shows a shoulder-type peak at lu A, which originates 
from the interband transition. 

Here, we examine the w-dependence of the Hall co- 
efficient Rh{'jj) = cfxyi'-^) / [c'xxi'-^))'^ and the Hall angle 

(lj) ~ axy{ijj) / (Jxx{^)- The real and imaginary part of 
Rh{'jj) is shown in Figs. [7](a) and (b), respectively. From 
these figures, we see that both Re (w) and Im Rh{oj) 
changes significantly at w ~ A for small 7. After reach- 
ing a peak at a; A, Re Rni^) changes its sign, whereas 
Im Rh{ijj) remains negative. However, this sharp peak 
at a; A can be easily suppressed as 7 increases, and Re 
Rh{'jj) for 7 = 0.1 and 0.5 remains negative. As for the 
Hall angle, its cj-dependence is shown in Fig. [51 When 7 
is small, Im 9h{uj) shows a peak at w ~ A, and it changes 
its sign for lu > A. 

Here, we briefiy discuss the w-dependence of Rh{(^) 
and 9h{(^) based on the simple Drude model. In the 

case of w < A, a^y{uj) - ct" —-^ ~ '^"^'^'-'^^ 



xy 



(0) and (Txxii^) 



1 — iuT 

where t = I/27. Then, Rh{uj) ~ (0)(1 - iur)'^ 



cx 



(47^ 



-UJ 



-Ai^uj, and Onito) ~ 0ij(O)(l — zwr) oc 27- 



lUJ. 



This analysis well explain the numerical results in Figs. 
[7]and[5]for w < A. 
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V. CALCULATION OF AC HALL 
CONDUCTIVITY WHEN 7 IS ENERGY 
DEPENDENT 

In the previous section, we calculated the AC AHC in 
the constant 7 approximation, assuming that the inelas- 
tic scattering due to local impurities are dominant. How- 
ever, in usual metals, inelastic scattering due to electron- 
electron interaction will be dominant, since the quasipar- 
ticle damping rate increases with lj: In a Fermi liquid, 
the imaginary part of self-energy ImE(ti;) is given by 



7(w) = -lmE(w) = /? (TrTf+u^ +7(0) 



(31) 



when T and u are not large, where /3 is a constant and T 
represents the temperature. 7(0) represents the damping 
rate at the Fermi level due to elastic scattering. Here, we 
study the AC AHE when the quasiparticle damping rate 
7(0;) is given by eq. (PT|) by putting 7(0) = 0.005. 

When the damping rate 7(0;) depends on uj, we cannot 
use eqs. (dH]) - (l^ni) . Therefore, we perform the numerical 



calculations in eq. (jl6p for axy{uj). To perform this, we 
decompose eq. (jl6p in the difference of integral interval 
as follows: 



(32) 



X G%,{e + {Gf,Ae - to/2) - Gj^e ~ ^/2)} , 

(33) 



-Lu/2 1 



h,a.,a.' ,13,(3' 

X {G%.{e + c./2)G^„,(e - u:/2) - {R ^ A)) . 

(34) 

Here, we remind the readers that only the terms 
{a',a,P\P) = {x,x,x,y),{x,x,y,x),{x,y,x,x) and 
(y, X, X, x) remain finite after fc-summation. Therefore, 
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A=0.417 




FIG. 9: w-dependence of Re axy{ijj) (/3 — 0.1, 0.5 and 1) when 
the dammping rate 7(1^-') is given by eq. pi|l . Here, Re a^cyii^) 
obtained for constant 7 (7 = 0.03, 0.05) are shown for com- 
parison. 



eqs. (|33p and ([34)1 are rewritten as follows: 



;/2 



uj + Z7(e + uj/2) - Z7(£ - uj/2) 

d«(e + w/2)d«(e-w/2) 
u) + i-f{e + tj/2) + i-f{e - uj/2) 



d'^{e + uj/2)d^{e~uj/2) 



(35) 



2X ^ 

TTLlI ^ 
k 



X ^ ti) X Im 



0/2 



d^{e + uj/2)d^{e-uj/2) 
+ {7(e + c^/2)-7(£-^/2)} 

xRe ' 



d^ie + Lu/2)dR{e -Lu/2) 



(36) 



We perform the fc-sumniations in eqs. psp and 
numerically, deviding the Brilliouin zone into 1000 x 
1000 meshes. As for the numerical e-integration, we per- 
form /y'^^^ de in eq. by setting X = 500. Since the 
domain of e-integration in eq. (|36p is about 100 times 
larger than that in eq. (j35p . e-integration is performed 
by dividing it into 50000 meshes for the former integra- 
tion, and 500 meshes for the latter integration. 

Now, we show numerical results. Figure [9] shows the lu- 
dependence of the AC Hall conductivity for (3 = 0.1, 0.5 
and 1, which corresponds to 7(A/2) = 0.0094,0.027 and 
0.049. The result of Re axyiut) for constant 7 (7 = 0.03 
and 0.05) are shown for comparison. We see that ob- 
tained Re axy{(^) using eq. ([?!]) is well reproduced by 
the constant 7 approximation for w ~ A, by putting 
7 = 7(A/2). To explain this fact analytically, we derive 
an analytical expression for Re a^yi^) when 7(0;) is co- 
dependent. At zero temperature, axyi^o) is mainly given 



by, 



axy{ijj) oc 



1 



1 



e + uj + fi-E+ + ij{e + Lu)e + ^i-E^^ - ij{s) 

1 1 
e + UJ + II - Ej^ + i^{£ + u!) e + fi- E^ - i^{£) 

(37) 

We note that the first term represents the interband tran- 
sition at w ~ E'^ — E^ , and the second term corresponds 
to ~uj ^ E^ — Ej^ . Hereafter, we neglect the second term 
since we study w > in the present study. By taking the 
pole of first term in eq. p7p . a^yi^) can be estimated as 



axyiuj) oc 



z r deJ20{E+ ~ ^i)e{^l ~ E^) 



'7r6{e + uj + fj.- E+) 



UJ — z7(e + uj) — i'y{e) 
TT5{e + fi-E^) 



-K + Ek + ^ + hie + Lu)+ Z7(e) 
Therefore, axy{uj) for uj ^ E^ — E^ is given as 



(38) 



(Jxv[Uj) OC I > -T- -T- ; -. 

^uJ-{E+-E-)+^^{E+-^i) + ^^{E--^i) 

(39) 

In the case of ~ A, the fc-summation in eq. (|39p is 
restricted around k'^. Then, we approximate as E^ ~ 
IJL + A/2 and E^ n - A/2 0, |^. As a result, we 
obtain the following extended-Drude (ED) expression for 
w - A: 



(Txyiuj) OC 



UJ - A + 2z7ed(w) ' 
where 7ed(w) is approximately given by 

7(A/2), 



(40) 



7ED(c^) = Tr[7(A/2)+7(-A/2)] 



(41) 



For example, when /3 = 1 and 7(0) = 0.005, the damping 
rate at w ~ A is estimated as 7(A/2) 0.049. We 
also verified that Re axy{uj) is mainly given by the Fermi 
surface term when the damping rate "f{uj) depends on uj. 

Here, we remark that eq. (j3ip is appropriate only for 
small UJ. According to the perturbation thoery with re- 
spect to Coulomb interaction, j{uj) will increase drasti- 
cally for > A due to the interband excitation. This 
fact will suppress axy{uj) for uj > 2 A further. There- 
fore, for a more reliable study, we have to calculate the 
w-dependence of 7(0;) microscopically. 

Finally, we comment on two important future prob- 
lems. The first one is the detailed study of the Coulomb 
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interaction effect on tlie AC AHE. According to the mi- 
croscopic Fermi liquid theory, the effect of Coulomb inter- 
action is exactly renormalized to the self-energy correc- 
tion and the CVC. As we have shown, the imaginary part 
of the self-energy tends to suppress the AC AHC. As dis- 
cussed in ref. [l6| , the renormalization factor due to the 

-1 



ex- 



real part of the self-energy, z = |^1 — 

actly cancels in the formula of the AC AHC given by eq. 
P^ . On the other hand, it is well-known fact that the 
CVC due to Coulomb interaction causes various anoma- 
lous transport phenomena in the vicinity of the magnetic 
quantum critical points (QCP) 62, 63, 64, 65, 6(|. This 
fact suggest that the CVC may cause novel temperature 
dependence of the AC AHC near the magnetic QCP. This 
is an important future problem. 

Another future problem is to perform a more reliable 
calculation on AC AHC based on a realistic tight-binding 
(TB) model. Recently, we have calculated SHC in Ad- 
and 5(i-transition metals based on the Naval Research 
Laboratory tight-binding (NRL-TB) model fUHl. This 
model enables us to construct nine- orbital {s+p + d) TB 
models for each transition metal [H, [lO]. In the future, 
we will study AC AHC based on this model to obtain 
more reliable results for AC AHC. 



VI. 



DISCUSSIONS 



A. Comments on Experiments 

In sections IIVI and |Vl we have discussed the 7- 
dependence of AC AHC. We found that the spiky peak 
at A exists when 7 is very small, whereas it is easily sup- 
pressed as 7 increases. The value of the damping rate 
7(0;) at cj ~ A/2 determines the characteristic behavior 
of the AC AHC. Here, we show that the spiky peak may 
vanish in d-clectron systems by using the experimental 
value of 7(ijj): Cerne et al. [61j reported the damping 
rate 7(0;) in Au and Cu. As for Cu (Au), the damp- 
ing rate 7(w) is obtained as ~ 600 (700) cm"-"^ when 
the frequency is ~ 900 (1000) cm^^. This means that 
7 « 600 K ^ 0.15 at w « 900 K - 0.25 in the present 
unit of energy. Since A > lOOOK in usual transition met- 
als [13, HI], the observed damping rate is large enough 
to suppress the spiky peak of the AC Hall conductivity. 
Furthermore, the damping rate in the strongly correlated 
systems will be larger than those in Au and Cu. There- 
fore, we conclude that the peak of a^yiuj) at w ~ A will 
be tiny or absent in usual transition metal ferromagnets. 

As shown in Figs. [4]and[5l the intrinsic AC AHC shows 
prominent deviation from the Drude-like behavior. On 
the other hand, the AC AHC due to the skew-scattering 
mechanism is expected to follow the Drude-like behav- 
ior, axy oc (7 — iuj)~^ for small ui. Therefore, AC AHC 
measurements will be quite useful to distinguish between 
the AHE due to intrinsic effect and that due to extrinsic 



effect, without necessity to introduce disorders. 

B. Summary of the Present Study 

In this paper, we studied the intrinsic AC AHE in 
transition metal ferromagnets based on {dxz, (ij,z)-orbital 
tight-binding model. We drived an analytical expression 
for the AC AHC that is valid for any quasiparticle damp- 
ing rate 7 without CVC, by performing the analytic con- 
tinuation carefully. We find that the intrinsic AC AHC 
does not follow the Drude-like behavior. When 7 is very 
small, AC AHC has a spiky peak at cj ^ A, which arises 
from the interband transition as explained in Fig. 3] (a). 
This behavior corresponds to the previously reported re- 
sults by Fang et al. 14]. When 7 is finite, however, the 
spiky peak is easily suppressed to be small or absent. In 
this case, the magnitude of AC AHC remains almost un- 
changed in the region < a; < A. We also calculated 
axx{i^), the Hall coefficient RH{i^), and the Hall angle 
9h{^jj)'- Rh{^) and 9h{uj) show a sharp peak at w = A 
for small 7, whereas this peak is easily suppressed as 7 
increases as shown in Figs. [7]and[51 

The overall behavior of the AC AHC axyiio) is repro- 
duced by the Fermi surface term (/) , as in the case with 
the DC AHC. The Fermi surface term strongly depends 
on 7, whereas the Berry curvature term (lib) has a weak 
dependence of 7. Although the relation axy{uj) « ct^^^((jj) 
holds in the present model when 7 is very small, the re- 
lation axy{oj) « axy{(^) is well satisfied for a wide range 
of 7. Therefore, prominent 7-dependence of AC AHC is 
well reproduced by the Fermi surface term. We stress 
that axy{i^) 7^ "'lii? (tiii even if 7 = +0 in general multi- 
orbital systems [l^, |43] . For a quantitative study of the 
intrinsic AHC, however, we have to calculate both the 
Fermi surface term and the Fermi sea terms on the same 
footing. 

Finally, we comment on the AC SHE. Recently Kon- 
tani et al. [i^ have studied the intrinsic SHE in d- 
electron systems. Therein, they have found that the 
present {dxz,dyz) tight-binding model shows huge SHE. 
In the present model, SHC ct^j, is given by {—h/e) times 



the AHC axy- o"™ 



;crxy, since the spin of the con- 



duction electron is conserved. Therefore, interesting in- 
dependence of AC AHC derived in the present study is 
also expected to be realized in AC SHC in various tran- 
sition metal complexes. 
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